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ABSTRACT 

The mutual and tripartite information between pairs and triples of disjoint regions in a 
quantum field theory are sensitive probes of the spread of correlations in an equilibrating 
system. We compute these quantities in strongly coupled two-dimensional conformal field 
theories with a gravity dual following the homogenous deposition of energy. The injected 
energy is modeled in anti-de Sitter space as an infalling shell, and the information shared 
by disjoint intervals is computed in terms of geodesic lengths in this background. For 
given widths and separation of the intervals, the mutual information typically starts at its 
vacuum value, then increases in time to reach a maximum, and then declines to the value at 
thermal equilibrium. A simple causality argument qualitatively explains this behavior. The 
tripartite information is generically non-zero and time-dependent throughout the process. 
This contrasts with (but does not contradict) the time-independent tripartite information 
one finds after a two-dimensional quantum quench in the limit of large time and distance 
scales compared to the initial inverse mass gap. 



1 Introduction 



Consider a quantum field theory whose state is described by the density matrix p. Mea- 
surements made within a spatial region A (with complement A) can all be described in 
terms of the reduced density matrix pA = tr^p. The entanglement entropy of A is defined 
as S{A) = — tr [p^ log Pa]- If the full system is in a pure state, Sa equals Sa and measures 
the amount of entanglement between A and A. Now consider two disjoint regions A and B. 
The amount of correlation between these regions (both classical and quantum) is measured 
by the mutual information 



While the entanglement entropy of a spatial region A has a UV divergence proportional 
to the area of the boundary of A, the mutual information is finite, making it an especially 
convenient quantity to study. The mutual information appears in bounds that limit how 
well an observer of A can predict events in B. Indeed I{A, B) > with equality if and 
only if A and B are uncorrelated, i.e. Paub = Pa® Pb- Mutual information can potentially 
provide a powerful description of how correlations evolve and spread in an out-of-equilibrium 
system. 

Progress in computing mutual information has been made for certain two-dimensional 
conformal field theories in special limits [Ij. In particular, Calabrese and Cardy have 
studied the evolution of entanglement entropy following a quantum quench j2|. They studied 
quenches in which a two-dimensional field theory with a gap suddenly becomes conformal 
(e.g., because a parameter in the Hamiltonian such as an external field is changed). The 
evolution of mutual information in this system can be computed in limits where the time, 
the interval sizes, and the interval separations are all large compared to the initial inverse 
mass gap. It was found that the entanglement entropy of an arbitrary number of intervals 
(and hence the mutual information defined above) exhibits a piecewise linear behavior in 
time. This behavior is consistent with a simple model based on causal propagation with 
the speed of light from an initial state with short-range correlations. 

Another interesting quantity to consider in this context is the tripartite information, 
which measures the degree of extensivity of the mutual information. This quantity is 
defined for three spatial regions A, B and C as 



h{A,B,C) = S{A) + S{B) + S{C)-S{AUB)-S{AUC)-S{BUC) + S{AUBUC) 



By definition, it is symmetric under permutations of its arguments. In a generic field theory, 
depending on the choice of the regions, it can be positive, negative or zero, meaning that 
the mutual information is subextensive, superextensive or extensive, respectively. It does 
not have the UV divergences present in the entanglement entropy. This is true even when 
the three regions share boundaries (in contrast to the mutual information which is cut-off 
dependent in this case and needs to be renormalized) . 

For equilibrium states of field theories that have a gravity dual, the AdS/CFT corre- 
spondence gives us a simple recipe for computing entanglement entropy in terms of the 
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areas of minimal surfaces in a dual spacetime (sl. The mutual and the tripartite informa- 



tion between regions can then be computed through (1.1 ) and (1.2). Applying this method, 
an interesting phase transition was observed, where the mutual information between a pair 
of disjoint intervals in two-dimensional field theory is of 0(1) for large separations, but of 
order the central charge for small separations f3|l4]. A recent work has also proven that 
for quantum field theories at equilibrium with a holographic dual, the mutual information 
is always extensive or superextensive, meaning that Is{A,B,C) < 0. 

The purpose of our paper is to understand how these results extend to a dynamical 
setting where energy is injected at t = and then proceeds to equilibrate. We model the 
injection of energy as a shell of null dust falling into anti-de Sitter (AdS) space. We then 
apply a proposal of [6] for computing entanglement entropy in a time-dependent AdS / CFT 
setting in terms of the areas of extremal surface in the dual spacetime. A key difference 
between our setting and that of [2] is that our theory is always conformal and hence starts 
with long range correlations, whereas their model is gapped at early times so that the 
correlations are essentially local [7|. 

We find that the mutual information between intervals in our setting starts at the 
vacuum value and ends at the thermal value, but usually passes through an intermediate 
phase where it is higher than either. The rise and fall are nearly linear in time. A causality 
argument (slightly modified from [2|[7]) explains qualitative features of our results. We find 
that the tripartite information is non-positive in all the cases that we study, matching the 
expectation based on the results of ^ who studied static configurations. In our dynamical 
setting the tripartite information changes with time during the process of thermalization in 
many parameter regimes. In contrast, the Calabrese and Cardy results for two-dimensional 
quantum quenches, where a mass gap is suddenly taken to zero, imply that is constant 
in their setting j2]. If the AdS/CFT entanglement entropy proposal of |6] is correct, the 
difference is probably due to the fact that our initial state has long distance correlations. 
In other words, the "mass gap" in our setup would be zero, so that we are studying the 
system on time and distance scales small compared to the would-be inverse initial mass 
gap. 

2 Mutual information 

Consider the 3-dimensional infalling shell geometry described by the Vaidya metric 

ds"^ = -[r^ - rjjOiv)] dv^ + 2dr dv + r'^dx^ , (2.3) 

where v labels ingoing null trajectories, is the step function and we have set the AdS 
radius equal to 1. In this form, the Vaidya metric describes a zero thickness shell composed 
of tensionless null dust (Fig. [l^). Outside the infalling shell (f > 0), the geometry is 
identical to that of an AdS black brane with Hawking temperature T = rjijli^^ while inside 
the shell [y < 0) it is the same as that of pure AdS. The change of coordinates 

t+^ln^ v>0 

2rH r+TH { 2 4 1 

t- i W < ^ ' ' 
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brings the two metrics to the standard form 



— r 



H 



)de + 



' H 



-r^dt^ + ^ + r^dx^ 



v>0 
V <0 



(2.5) 



On the boundary, at r = oo, the coordinates v and t coincide. 



singularity 




Figure 1: (a) The causal structure of the Vaidya spacetime shown in the Poincare patch 
of AdS space. The asymptotic boundary (vertical line on the right hand side) is planar, 
and the null lines on the left hand side of the diagram represent the Poincare horizon, 
(b) Connected (in blue) and disconnected (in red dashed) locally minimal surfaces for the 
boundary region AU B in AdSa. 



We consider the time evolution of the mutual information between disjoint intervals 



of length i separated by a distance d on the boundary of the Vaidya spacetime (2.3). In 
three bulk dimensions the entanglement entropy of a connected region A on the spacetime 
boundary is proposed to be computed holographically by the length of the bulk geodesic 
that connects the endpoints of A |6|. For disjoint intervals A and B there are multiple 
candidates for the geodesies that contribute to S{A U B). In [2], it was shown that for 
equilibrium configurations, S{A U B) is given by the minimum of the lengths of two sets of 
geodesies (shown as solid and dashed lines in Fig. [i]d). We will assume that this prescription 
extends to the dynamical Vaidya background as a rule saying that S{A U B) is determined 
by the length of the shortest collection of geodesies connecting the endpoints of A and B^ 
The geodesies connecting two equal-time boundary endpoints (to, 0) and (to, ^) have been 
studied in (Tjlsl. These results have been extended to non-equal time spacelike geodesies 



in 19 . Since the geodesic length C diverges due to contributions near the AdS boundary. 



-•^ There is also a topological condition Q that the geodesies should be continuously deformable to the 
AdS boundary. This condition will not come into play for us. 



3 



it is more convenient to work with a regularized length computed by introducing a large 
radius cutoff in AdS at r = tq and defining bC = C — 1 ln(2ro). This removes the divergent 
part of the length in the vacuum state (pure AdS). For to ^ (before the shell falls in) the 
geodesies lie entirely in AdS at fixed time t = to. Their renormalized length is 



-^vacuum (^) 2 In 



(2.6) 



For < to < i/2, the geodesies cross the infalling shell and their renormalized length is 

sinh(ri^to 



^£sheii(^,to) = 21n 



rHs{£,to) 



(2.7) 



where s(£, to) ^ [0, 1] is parametrically defined by 



1 



2c /2(l + c)p2 + 2sp-c 
Jp^ V2(l + c)p2-2sp-cy_ 



(2.^ 



with c = Vl — and 



2p = coth(ri^to) + \ coth {rnt 



2c 



(2.9) 



For to > i/2, the geodesies lie on the t = to surface in the black brane geometry and 



thermal 



sinh^ 
2 In ^ 



(2.10) 



In terms of these geodesic lengths our formula for the mutual information between the 
intervals A and B is 



/(A 5) 



1 



AG 



N 



{25C{i) - Min [25£(£), 5£(2£ + rf) + bC{d)\} . 



(2.11) 



Observe that in principle we should also consider geodesies of spatial separation £ + 
which connect the two left (right) end-points of the two intervals A and B. We should thus 



evaluate in (2.11 ) 



Min [25£(£), bL(2t + rf) + bL{d), 2bL{i + d)\ . 
However, for any fixed boundary time to: 

-^vacuum (^ ~l~ C^) ^ -^vacuum (^) ; 
^'Cthermal(^ + C^) > '^'^thermal (^) 5 
5£shell(^ + to) > 5£shell(^, to) , 



(2.12) 



(2.13) 



since, in each case, bL is monotonically increasing as a function of the spatial boundary 



separation, as is seen from (2.6), (2.10) and numerically from (2.7). Therefore, when all the 



possible competing geodesies are of the same type (in terms of crossing the shell or not), it 
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Fig lire 2: (ithemial cUld (^vacuum &s a function of £ for r^^ = 1. The mutual information surely 
vanishes for to < and to ^ (2^ + (i)/2 for d and i in the white region, for to > (2^ + d)/2 
in the green shaded region, while it is generically everywhere non-zero in the blue dotted 
region. 



is enough to consider (2.11 ). There are however time intervals in which the property (2.13) 
is not enough to exclude the case of intersecting geodesies of length 6C{i + d) from (2.12). 
For example, when £ > d and to G [£/2, {£ + d)/2), one has 6C{i + d) = 5£sheii(^ + c?, to) 
while 5C{£) = 5Ahermai(^)- In these cases, however, it follows from ( 2.7[ ) and (2.10) that 



(5£shell(^ + d,to) > (5£thermal(^) 



(2.14) 



in the relevant time interval, and it is therefore always enough to consider (2.11). 

In the Vaidya spacetime (2.3), for to < 0, the mutual information of A and B coincides 
with the vacuum AdS result of [4]: 



I{A,B) 



vacuum 



fin 



d(2£+d) 



d> d^ 
d<d^ 



(2.15) 



vacuum 



where we have used (2.6), the central charge c = 3/{2G]\f) and (ivacuum(^) = {V2 — 1)^- For 
^0 ^ (2^ + d)/2, all the geodesies lie entirely in a black brane background and the mutual 



information is given by 10 
1 



I {A, -B) thermal 



AG 



N 



{25£thcrmal(^) - Miu [25£thcrmal (^) , 5Ahermal(2^ + d) + 5£thermal(c?)] } 



In 



0, 

sinh^(rj:j£/2) 

sinh(rH(2^+d)/2) smh(rHrf/2) 



d ^ C?thermal 
d ^ '^thermal 



(2.16) 



using (2.10), where 



(^thermal (^) = In 



1 - e'"^ + e^'"^ + 7(1 - e'-H^ + e^'H^)"^ - e^'^i 



2i. 



(2.17) 
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Figure 3: Rescaled mutual information / = AG^I as a function of boundary time to for 
d = 0.4 (left), d = 2 (center) and d = 4 (right), and = 1. The various curves correspond 
to different values of i which increases from the bottom up. The left panel shows i = 
0.2, 0.4, . . . , 2.0, while the center panel and the right one show £ = 1, 2, . . . , 10 (some of the 
curves in the three panels are not visible because everywhere vanishing). 



As the intervals size i ^ oo, (ithermai(^) d^o = 

From these results, we see that for a fixed interval size i when ^thermal < d < livacuum the 
mutual information must vanish for to > {2i + d)/2. The reason is that at these late times, 
all the geodesies that can potentially contribute to the mutual information (see Fig. [T)d) will 
lie in the black brane background. When d > c?vacuum it is zero also for to < because at 
early times all the relevant geodesies lie in the empty AdS background. When d < ^thermal 
the mutual information can be non-zero for all times (Fig. |2]). 

In the Vaidya background the mutual information should interpolate between the vac- 
uum and thermal results. The interpolation occurs in the period < to < {2i + d)/2, when 
one or more geodesies intersect the shell. Crossing the shell has the effect of lowering (rais- 
ing) the geodesic length for a fixed boundary separation, with respect to the black brane 
(AdS) background. When Max(£/2, d/2) < to < {2i+d)/2 only the outermost geodesic with 
boundary separation 2i + d crosses the shell. When Mm{£/2,d/2) < to < Max{i/2,d/2) 
geodesies with boundary separation 2£ + d and Max{£,d) cross it, and so on. The mutual 
information in this range of times can be computed using the appropriate geodesic lengths 
SCsheii in (2.7) and 5£thermai in (2.10). The mutual information computed in this way starts 
at the vacuum value at early times and ends at the thermal value at late times. Remark- 
ably, it increases sharply to a peak at intermediate times (Fig. |3]). These results are easily 
extended to a situation where the two intervals have different lengths ii and £2- The logic 
is identical, so we merely plot the results in Fig. |4j 



3 A qualitative explanation for the mutual information? 

Why does the mutual information increase and then decrease again in our setup? For 
understanding the entanglement entropy of a single interval in the same setting it was 
useful to compare the AdS results with exact computations for quantum quenches [2||7j. 
The entanglement entropy for multiple intervals after a quench was also computed by 
Calabrese and Cardy |2| for times and spatial scales large compared to the inverse mass 
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Figure 4: Rescaled mutual information / = AG^I as a function of boundary time to for 
d = 0.4, £i = 2 (left), d = 2, £i = 8 (center) and d = 4, £i = 8 (right), and th = 1- 
The various curves correspond to different values of £2 which increases from the bottom 
up. The left panel shows £2 = 0.4, 0.8, . . . , 4.0 while the center and the right panels have 
4 = 2, 4,..., 20. 



gap of the early time state. From their results we can compute the mutual information 
following (1.1). Since the Calabrese-Cardy theory starts, unlike ours, with a gapped initial 



state and a type of energy deposition with short range correlations, their results need not 
agree with ours ||7]|. However, the explicit result nevertheless has a sharp feature in the 
middle. 

A simple toy model developed by [2] gives a physical interpretation for this principal 
feature of the mutual information in the quantum quench. The idea was that before the 
quench there is a mass gap and a finite correlation length. Thus, after the quench only 
excitations originating at almost the same point will be entangled. At later times, such 
a pair of correlated particles will only contribute to the mutual information between two 
intervals, if each interval contains one of the two particles. As shown in Fig. |5] this naive 
model would lead to a linear rise and subsequent decline in the mutual information. An 
extension of this toy model to the Vaidya setup reproduced the entanglement entropy of a 
single interval by including power law long range correlations in the initial state [Tj. 





(a) (b) 

Figure 5: (a) A simple causality derived picture for two intervals after a quench. After the 
quench at t = signals propagate outwards at the speed of light in both directions. These 
will give correlations which contribute to the mutual information at t = to if one signal is 
in each interval, which occurs only for signals originating in the shaded pink diamond, (b) 
The analogous picture for two intervals of different length. 
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Figure 6: Rescaled tripartite information J3 = AG^h for three adjacent intervals of length 
i, d and £ as a function of boundary time to, for d = 0.4 (left), d = 2 (center) and d = 4 
(right), and th = 1. The various curves correspond to different values of ^. The left panel 
shows ^ = 0.2, ... 1, in continuous lines with i decreasing from the bottom up, while the 
curves in dot-dashed have i = 1.2, . . . , 2.0 increasing from the bottom up. The center and 
right panels have £ = 1, 2, . . . , 10 decreasing from right to left. 

To try to reproduce our results more accurately, further modifications of the toy model 
would be needed. Consider, for example, the non-zero final state mutual information that 
is present in some cases. This can occur in a naive particle picture if the initial energy 
deposition had a finite spatial correlation length, and we consider particle pairs that travel 
in the same direction as well as in opposite directions. In this setup, two intervals separated 
by a distance d will always be populated at late times with particle pairs that move in the 
same direction and were produced at a separation slightly larger than d. If the initial 
correlation length was bigger than d a final state mutual information can result. 



4 Tripartite information 



Now consider the time evolution of the tripartite information (1.2) for three intervals on 



the boundary of the dynamical Vaidya background (2.3). From (1.2), in addition to the 
computations that we have already done, we need to compute S{A U B U C). Following 
the previous sections, we will assume that this is determined by the length of the shortest 
collection of geodesies connecting the endpoints of A, B and C. 

We first consider the case of three adjacent intervals A, B and C, of length £, d and 
£ respectively, on the boundary of the Vaidya spacetime. This corresponds to the setting 
considered in Sec. |2] for the mutual information, now promoting the region of length d, 
which separates the two intervals of length £, to a third interval. 



We plot in Fig. [6|the results for the tripartite information (1.2) for different values of 
£ and c? as a function of the boundary time to- For large values of d (center and right 
panels of Fig. [6]) the tripartite information starts at the vacuum value and ends at the 
(larger) thermal value passing through an intermediate phase where it is more negative 
than either. For small d (left panel of Fig. [6]) and large enough £, the behavior is different. 
In this case, the vacuum value is less negative than the thermal one and the tripartite 
information just decreases non-linearly. Incidentally, in this parameter range the thermal 
mutual information is non-vanishing (see Fig. Isl and values of ^ and d in the blue dotted 
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region of Fig. |2j) . 

For the value of d in the center panel of Fig. |6} for some values of the tripartite infor- 
mation goes through an intermediate phase where it is constant. This can be understood 
looking at the first line of (1.2). In this phase, all the contributions are constant except 



for S{A U C) and S{A U 5 U C), whose time-dependent contributions cancel each other in 



( |L2| ). To see this, first note that S{A), S{B), S{C), S{AUB) and S{BUC) all correspond 
to geodesies extending in the thermal background and have a constant value. Meanwhile 
S{A U B U C) is given by the geodesic that ends on the two endpoints oi AU B U C, which 
crosses the shell. Finally, S{A U C) is determined by the two geodesies enclosing a con- 
nected region in the bulk - the outer geodesic of this pair is the same one that appears in 
S{AU B UC). Therefore the two contributions that cross the shell and are time-dependent 
cancel each other between S{A U B U C) and S{A U C), leading to a constant result for 
the tripartite information. For large enough i/d, this constant value is approximately i- 



independent, as can be seen from (2.10). The same type of argument applies to late times. 



when all the geodesies extend in the thermal background. 

As shown in Fig. |6} we find a non-positive tripartite information at all times. Even 
though we are dealing with a dynamical setup, this behavior matches [s], where it has been 
shown that in holographic theories at equilibrium the mutual information is extensive or 
superextensive. 

While the case of adjacent intervals is a direct extension of the computation of the 
mutual information (Sec. [2]), the tripartite information for three disjoint intervals is more 
subtle. However, using the fact that only non-intersecting geodesies are local minimal 



surfaces for two disjoint intervals, as was the case in (2.11), the collection of geodesies that 



potentially contribute to the tripartite information of AU BUC in the Vaidya background 
reduces to those shown in Fig. [7} 



Figure 7: A schematic representation of locally minimal surfaces for the boundary region 
y4 U -B U C in AdSa. The top right figure with A and C interchanged also contributes. 
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The time evolution of the tripartite information for three disjoint intervals of length £, 
each separated by a distance d from adjacent ones, is plotted in Fig. [8] for various values 
of d and i. It mainly shows the same features as the tripartite information of adjacent 




Figure 8: Rescaled tripartite information J3 = AG^I^ as a function of boundary time 
for d = 0.4 (left), d = 2 (center) and d = 4 (right), and th = 1- The various curves 
correspond to different values of i which decreases from right to left. The left panel shows 
i = 0.4, 0.8, . . . , 4.0, the center and right panels have i = 1, 2, . . . , 10 (although only £ = 
3, . . . , 10 and £ = 6 . . . , 10 are not everywhere vanishing in the center and right panel, 
respectively) . 



intervals. 

We have found that the tripartite information computed in the Vaidya geometry is 
generically non-zero. Now recall that Calabrese and Cardy found that the entanglement 
entropy of N disjoint intervals (M2j-i, ^2j), with 1 < j < and Uk < Uk+i, after a quantum 
quench [2] evolves in time according to the formula 

2N 

^ ~ 5oo + ^ V (-l)'=-'-^Max[Mfc -t,ui + t]. (4.18) 

Here c is the central charge and tq is a parameter of the same order as the inverse mass gap 
(refer to [2] for more details). This formula holds only in the limit where the time t and 



all separations \uk — ui\ are much larger than the initial inverse mass gap. Using (4.18), it 
follows that the time-dependent part of the tripartite information vanishes in this regime 
because each contribution associated to an extremum uj appears an equal number of times 



with a minus and a plus sign in (1.2). In contrast, the tripartite information for the Vaidya 



background (Fig. ^ is generally time- dependent and is only constant (and vanishing) if 



I {A, B), I {A, C) and I {A, U C) in are all identically zero. 

The results of Calabrese and Cardy need not agree with ours since they start with a 
gapped initial state which lacks long-range correlations. Nevertheless it is interesting that 
important qualitative features of the time evolution of the mutual information matched 
between their model and ours, while the behavior of the tripartite information is sharply 
different. In any case, the disagreement shows clearly that the simple causality argument 
of Sec. [3} where interactions were not taken into account, is not enough to capture the 
qualitative features of our strongly coupled two-dimensional dynamics. 
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